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We propose a heterogeneous reaction-diffusion model for the hydration and setting of ce- 
ment. The model is based on diffusional ion transport and on cement specific chemical dissolu- 
tion/precipitation reactions under spatial heterogeneous solid/liquid conditions. We simulate 
the spatial and temporal evolution of precipitated micro structures starting from initial random 
configurations of anhydrous cement particles. Though the simulations have been performed for 
two dimensional systems, we are able to reproduce qualitatively basic features of the cement 
hydration problem. The proposed model is also applicable to general water/mineral systems. 
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I. INTRODUCTION 

In the present paper we propose a heterogeneous 
reaction-diffusion model for the hydration and setting 
of cement. The proposed model is based on the ex- 
perimental observation that cement hydration can be 
described by a dissolution/precipitation mechanism 

01- 

The elementary aspects of the cement hydra- 
tion/setting process on a mesoscopic length scale can 
be characterized as follows. At the initial stage ce- 
ment particles or powder (in our case tricalcium sili- 
cate, C3S = Ca^SiO^) are mixed well with the solvent 
(water). Rapidly after mixing the dissolution reac- 
tion of the cement particles starts. Its principal re- 
action products are ions which arc mobile and may 
diffuse into the bulk of the solvent (in our case Ca 2+ 
, OH~ and H 2 Si0 2 ~ ions). At this stage the ion 
concentrations in the bulk of the solvent are very low 
and, as a consequence, one finds strong ion fluxes from 
the dissolution front away into the solvent. However, 
for a given temperature and pressure the ion concen- 
trations can not take arbitrary high values. Rather, 
the ion concentrations are bounded by finite solubil- 
ity products above which solid phases start to pre- 
cipitate from the solution. There are two associ- 
ated precipitation reactions: a) the precipitation of 
calcium hydro-silicate sometimes referred as 'cement 
gel' (d. 5 SH 2 . 5 = (GoO)i.b (Si0 2 ) (H 2 0) 2 . 5 ) and b) 
the precipitation of calcium hydroxide or 'Portlandite' 
(CH = Ca(OH) 2 ) |§. While the growth of the ce- 
ment gel is the basis for the whole cement binding 
process, the growth of Portlandite mainly happens in 
order to compensate for the accumulation of Ca 2+ and 
OH~ ions in solution (see below). 

The process of cement dissolution, ion transport, 
and cement gel/Portlandite precipitation is usually 
referred as 'cement hydration'. This process is to a 
high extent heterogeneous in the sense, that starting 
from initially random cement particle positions disso- 
lution and precipitation reactions change the physical 
(boundary) conditions for the ion transports them- 
selves. This is due to the fact that the diffusion within 
the solid phases can be neglected relative to the sol- 



vent one. 

Coupled chemical dissolution/precipitation reac- 
tions are well known since a long time in geology and 
geochemistry. As such kind of processes are funda- 
mental for water/mineral systems some effort has been 
undertaken to modellize such systems |||7j]. The ap- 
proaches differ methodologically and focus on differ- 
ent physical and chemical aspects. Ref. [| considers 
the ion transport problem for a single, solubility con- 
trolled, dissolution/precipitation reaction employing 
a one dimensional cellular automaton approach. A 
chemically more detailed approach of water/mineral 
interactions includes considerations on the influence 
of nucleation, electro chemistry and temperature on 
chemical reaction rates, however, without considering 
a transport equation (7). 

Recently, a relatively simple stochastic cellular au- 
tomaton model for cement hydration/setting has been 
proposed || . Therein it has been tacitly assumed that 
ion transport is not relevant for the hydration process 
above a certain mesoscopic length scale ('pixel size'). 
Furthermore, in this model hydrate 'pixels' perform 
a random walk until they touch another solid 'pixel', 
where they then stick Jj| . In contrast to this we posi- 
tion our deterministic model into the 'opposite' direc- 
tion: mass transport happens due to diffusion of ions, 
and the hydrates are regarded as immobile. 

In Sec. H we give a detailed description of the em- 
ploy ed m odel. We consider heter ogene ity aspects in 
Sec. II A| , chemical aspects in Sec. [IB and transport 



aspects in Sec. II C. Sec. 



Ill contains first results. 



We present some images of calculated cement micro 
structures, which comprehensively demonstrate basic 
features and capabilities of the present approach. We 
consider parametric plots for the mean ion concentra- 
tions, known in cement literature as 'kinetic path ap- 
proach' and compare our findings with experimental 
results (9J. Following this we investigate the varia- 
tion of the maximum average silica concentration in 
solution for various chemical reaction rate constants. 
Finally we present curves for the hydration advance- 
ment in time for two different sets of reaction rate 
constants. In Sec. IV we summarize and give some 



ideas about future work. 



II. THE MODEL 

In the following we describe the proposed reaction 
diffusion model for the hydration and setting of ce- 
ment in more detail. 

We will first discuss how to quantitatively describe 
the heterogeneities and outline some general features 
of the employed model. Following this we will sum- 
marize the dissolution/precipitation reactions and the 
employed reaction rate laws. Finally some physical as- 
pects of the ion transport in the solvent are discussed. 
The later one represents the physical coupling between 
the dissolution and precipitation reactions. 



A. Heterogeneity aspects 

We consider a discrete reaction-diffusion model in 
space and time. The physicochemical system can be 
regarded as being composed of 'sufficient' small vol- 
ume elements V = Ax 3 with Ax typically between 
10 _6 m and 10~ 4 rn. The initial configuration is deter- 
mined from a digitized micro graph image as follows: 
the micro graph is immediately taken after mixing at 
time to j one can assume to find only the solid phase 
C3S and water. By measuring the occupied volumes 

of the phase a, (a = aq, C3S), in each cell i us- 
ing an image processing system, one determines the 

initial distribution of mole numbers n\ (to) accord- 
ing to nf\t Q ) = y. (Q) (t )/v (Q) . The ijW are known 
molecular volumes for room temperature, see Table 
I. In general, such a micro graph shows cement par- 
ticles of different sizes immersed in water. Because 
each cell volume is completely filled with solid and 
solvent phases, one can calculate the solvent volume 
of cell i at any time step from the corresponding solid 
volume (s), 



v l {aq \t) = v- £ 



4 a) (t)v^ a \ 



(2.1) 



This is a very equation, because ion concentrations 
are usually calculated with respect to the actual sol- 
vent volume Vj ^(t). The physical constraints of 



Eq.(2.1) are that the chemical reactions must be suf- 
ficiently slow in comparison to the solvent flow. Fur- 
thermore, it is tacitly assumed that the whole system 
is connected to an external water reservoir in such a 
way that the solid phases do not hinder the flow. 

In a given volume element the dissolu- 
tion/precipitation reactions will, in general, not in- 
stantaneously go to completion. The solvent volume 
in each cell will rather increase/decrease continuously, 
defining a solvent distribution field, 



S a <i) 



(t) = v^\t)/v. 



(2.2) 



The solid phase volume fractions e\ (t) are defined 
analogously. 

We believe in fact that the volume fractions e[ a \t), 
reflecting the actual solvent/solid distribution in the 
system, are sufficient to characterize most of the as- 
pects of the hydration process on a macroscopic length 



scale. 

Though we do not have information about these dis- 
tributions on a length scale smaller than Aa; we will 
assume that the reactants are homogeneously dis- 
tributed in each volume element. The volume frac- 
tions may then be interpreted as probabilities to find 
at a random position within cell i the reactant a. 



B. Chemical aspects 

In the following we will summarize the considered 
dissolution/precipitation reactions and the employed 
reaction rate laws. During the course of the simulation 
it may happen that the reaction flux, as calculated 
from the kinetic equations, does exceed the available 
amounts of chemical reactants in a given cell. We 
therefore check at each time step and for each chemical 
reaction if there exists a limiting reactant. If this hap- 
pens, we define the reaction flux through the amount 
of the available limiting reactant. 



1. C3S dissolution 

The C3S dissolution reaction is a spontaneous and 
exothermic surface reaction which happens at the 
CsS'-solvent interface(s) ||. It can be considered as 
irreversible, 



CsS( 3 ) + 3H 2 0(e) 



AOH, 



(aq) 



H2Si0^ aq y 

(2.3) 



Here k± denotes an appropriate surface reaction con- 
stant in units of molm s . 

Chemical reactions must be considered as local, i.e., 
each volume element is acting as a small and indepen- 
dent chemical reactor as long as no transport occurs. 
However, for a pronounced dissolution reaction in a 
given cell one has also to consider the possible reac- 
tion amounts originating from the interfaces with its 
neighboring cells. This can be done by allowing ce- 
ment dissolution in cell i through the electro chemi- 
calsolvent of that cell and through a fraction 7 of the 
solvent of the neighboring cells j, 



i a) dfci (t) = "1°° fc i At Ax2 (pa + 7 J2 (P*i -Pa)) 



,(«) 



j=nn(i) 



(2.4) 



Here v{ ' denote the stoichiometric numbers of 
species a in reaction (|]^), see Table I. The changes in 
mole numbers of species a in cell i in a time intervall 
At due to this reaction are denoted by v[ a ^ d£i t i (t). 
The sum on the right hand has to be taken over all 
next nearest neighboring cells of cell i, as indicated by 



j = nn(i). The first term p it = ef^e^ in Eq.fl2.4D 
describes the dissolution within cell i (7 = 0) due to a 
'typical' reaction interface. It can be understood as a 
probability to find the two reactants, C3S and water, 
in contact at an arbitrary chosen point within cell i. 
The 'chemically active interface' between cell i and j 



(C 3 S)(aq) 



The constant 



is given by pijAx 2 = Ax 2 e. 
7 controls the degree of dissolution due to this inter- 
face^). We have used a value 7 = 1/8 in our simu- 
lations, yielding comparable reaction rates from inner 
and outer cell surfaces. The determination of the re- 
action interfaces in terms of cement and water volume 
fractions is similar to the degree of surface coverage 
used in Langmuir absorption theory [ flO| ]. However, 
the absorption of water on the cement's surface and 
the desorption of ions from this surface are not nec- 
essarily the rate d eter mining steps. In fact, the dis- 
solution reaction (2.2) involves various physical and 
subprocesses [|| which will not be considered here. 



2. CSH dissolution/precipitation 

The precipitation of C1.5SH2.5 is an endothermic 
reaction which can be considered as reversible 01 , 



The stoichiometric numbers again are denoted by 
, see Table I, and the change in mole numbers of 



,(«) 



species a due to this reaction is d^i^{t) 



3. CH dissolution/precipitation 

The precipitation of calcium hydroxide accompa- 
nies the CSH precipitation. This is because the CSH 
precipitation does not consume the ions in the same 
proportions as they are released due to C3S disso- 
lution. The non-reacted Ca 2+ and OH" ions soon 
begin to accumulate in solution, until the solubility of 
CH is exceeded. However, the precipitation of CH 
and CSH are, in general, not simultaneous because 
the corresponding solubilities are very different, comp. 
Table I. The dissolution/precipitation of CH can be 
considered as reversible, 



C1.5SH2. 



5(s) 



Mo 4 2 ( - 9) 



CH, 



OH. 



(aq) 



(2.5) 



Ca 2+ ,+20H7 

(aq) (aq)' 



(2.7) 



having a dissolution constant k% {mol m 2 s 1 ) . 

The corresponding solubility product for CH is de- 
fined as Sl qm = S u \^ ui = (Ca 2 +),(Oi/-)f| e «™ » 
3.3 x 10~ 6 mol 3 liter" 3 which is about three orders of 



magnitude larger than S 2 , see Table I, 



4 a) d&,3 (t) = 4 a) h At Ax 2 e ™ (1 - ^ 



). (2. 



In all other respects the CH reaction is treated 
analogously to the CSH dissolution/precipitation re- 
action. 



4- Chemical shrinkage 



One can write down from Eqs. (|2.3|) , (2.5) and (2.7) 
the net reaction for the cement hydration process, 



having a forward (dissolution) rate constant ki (in 
units of mol m" 2 s" 1 ). 

We note, that the CSH precipitation, i.e., the back- 
ward reaction in Eq. (|2.5| ) is the main reaction in this 
balance equation. The reactants and products are 
considered to have a fixed stoichiometry. However, it 
is known from experimental data that the stoichiome- 
try of C x S y H z can be rather variable ('solid solution') 
We will not include these complications into the 
model. Instead we consider a fixed and typical cal- 
cium/silica ratio equal to 1.5. 

In equilibrium the C1.5SH2.5 dissolu- 
tion/precipitation reactions are controlled by an 
empirical value of solubility, giving the maximum 
amount of solid one can dissolve in aqueous solu- 
tion at a given temperature (room temperature) 
and pressure (normal pressure). Employing this 
empirical solubility constant one can directly de- 
termine the equilibrium solubility product, S^ 9 " 1 = 
S 2 i\ equi = {Ca 2 +)l- 5 (OH-)i(H2Si0 2 -)i\ e i ui rj 1.2 x 
1CT 9 mol 3 5 liter" 3 - 5 , comp. Table I. The square C 3 S {s) + 4ff 2 (£) C 1 . 5 SH 2 . 5 ( S ) + l.5CH {s) . (2.9) 
brackets denote the ion concentrations with respect 

to the available solvent, i.e., (Ca 2+ )i = n\ Ca ' /V^' 
etc.. When in a given cell the ion product is larger 
than the solubility product, locally precipitation hap- 
pens. If it is lower, C1.5SH2.5 becomes dissolved. The 
CSH dissolution is only marginal and for this reason 
we will employ for the rate of dissolution a simpler 
expression than for the C3S dissolution reaction. 

We assume that the dissolution reaction in cell i is 
proportional to the typical reaction interface Ax 2 Ej ' 
with the dissolution constant k% as the constant of pro- 
portionality. The precipitation reaction is assumed to 
be proportional to the reaction interface and propor- 
tional to the ion product S 2 i, as defined above. The 
rate constant of the precipitation reaction can be elim- 
inated, by use of the equilibrium condition, d^i,2 = 0. 
We have employed the following rate equation, 



d4i, 2 (t) = vf k 2 AtAx 2 e[ aq) (l 



>2i 



5equi , 
2 



(2.6) 



Calculating for the above equation the occupied vol- 
umes for one mole one finds for the left side approxi- 
mately 144 cm 3 and for the hydrate products 130 cm 3 , 
comp. Table I. Hence, the reaction products occupy 
a volume around 10% smaller than the reactants (in- 
cluding solvent). This effect, which is typical for hy- 
draulic binders, is termed 'chemical shrinkage' or Le 
Chatelier effect. It is responsible for the transition 
of cement suspensions or pastes having a finite vis- 
cosity and zero elastic moduli to rigid hydrated ce- 
ment having an infinite viscosity and finite elastic 
moduli. Wherever topologically possible water flow 
will try to compensate for the loss of volume. How- 
ever, in the course of an experimental cement setting 
and hardening (rigidification) process available sol- 
vent flow paths may vanish and volume loss becomes 
partially compensated building up mechanical defor- 
mations (stresses) within the solid phase agglomer- 
ate. In turn these stresses influence the reaction rates 



(which in general are pressure dependent) in a com- 
plex manner The appearence of voids (pores) on 
the submicron scale is experimentally also well estab- 
lished §. 

We will not modellize such complex problems as me- 
chanical stresses, solvent cavitation and solvent flow 
due to chemical shrinkage in our present approach. 

Instead we will assume throughout all calculations 
that the solvent is able to balance the chemical shrink- 
age for all tim e ste ps and for all volume elements ac- 
cording to Eq.( |2.l| ). This assumption is equivalent to 
the introduction of local source terms for water, i.e., 
the hydrating system is connected to an external wa- 
ter reservoir in an appropriate way. 



We would like to point out that Eq. ( 2.9 ) is an overall 
reaction which per definition only holds for an isolated 
system. On the other hand one cannot treat the vol- 
ume elements as chemically isolated, simply because 
of the ion fluxes going through the elements. This and 
the different values for the CSH and CH solubilities 
imply in general locally non-congruent precipitation 
reactions. As a consequence one cannot simply 're- 
place' one dissolved unit volume of C3S by 1.7 unit 
volumes CSH and 0.6 unit volumes CH as has been 
proposed in [||. 



The general formulation of the ion transport prob- 
lem in cement chemistry involves various subproblems 
which we will not consider at this stage of modeliza- 
tion, as there are the transport due to 



1. Heat conductance 



The cement dissolution reaction is strongly exother- 
mic. Analogous to the concentration gradients one 
finds near the various dissolution fronts strong tem- 
perature gradients, leading to non isothermal solvent 
flow fields, with its complications. The heat redistri- 
bution is probably also of direct importance for the 
dissolution/precipitation reactions, because their re- 
action constants as well as the corresponding solubil- 
ity products usually depend very sensitively on the 
temperature. We will neglect this possible effects in 
our approach assuming an overall constant room tem- 
perature, Tq = 298 Kelvin. 



2. Electrostatic interactions 



C. Transport aspects 

It is known from experiments that ion transport due 
to diffusion is fundamental in cement chemistry, be- 
cause of very strong ion concentration gradients close 
to the dissolving CsS'-solvent interface (s). In general 
one should consider the convective diffusion equation, 
however, this would imply to solve the hydrodynamic 
equations. In the present paper we consider the usual 
diffusion equation with local, appropriately defined 
transport coefficients. It is in the context of the inves- 
tigated problem useful to consider only the ion diffu- 
sion within the solvent. For convenience we document 
in Table I the diffusion constants as calculated 

from the electric mobilities at room temperature and 
normal pressure for infinite dilution. The solids and 
the solvent's diffusion constants are set to zero JjT[ . 
This allows a relatively compact notation for the equa- 
tions of continuity, 



» _ 



At 
Ax 2 



(«)„(«) _ nM/.teh 



j=nn(i) 



(2.10) 



fc=i 



The second term on the right side of Eq.( 2.10 ) rep- 
resents the already kno wn so urce term due to the 
chemical reactions Eqs.(2.4), ( |2.6| ) and (2.8). The 



are the transport coefficients, that depend on 
the solvent volume fractions (porosities) of the cells in- 
volved in the diffusion process, — D^ef q 'e^' q \ 
The ion concentrations are taken with respect to the 



actual solvent volume, i.e., c-"' 1 = n y "' /V^""" . We 
have employed periodic boundary conditions in all cal- 
culations. 



Because one is dealing with the transport of elec- 
trically charged particles, one must consider, in gen- 
eral, Poisson's equation for the electrostatic potential 
in order to determine the resulting ion mobilities due 
to their electric field. We will neglect in our simula- 
tions all electrostatic interactions. Furthermore, all 
chemical rate calculations are done in terms of molar 
concentrations and not in terms of 'chemical activities' 
(which is equivalent in the limit of infinite dilution). 
Hence we will treat the ions in solution as uncharged 
particles throughout all calculations. 



3. Solvent flow 



This is possibly of importance for the overall hydra- 
tion process, because the solvent has to move during 
the dissolution/precipitation reactions accordingly, in 
order to compensate for the loss/gain of solid volume, 
creating advective ion fluxes. For an infinite, plane 
cement-water interface one might argue that the sol- 
vent flow is always perpendicular to the interface. Un- 
der this particular conditions it can be shown by scal- 
ing arguments that advective ion transport is unim- 
portant relative to the diffusion transport, because of 
the relatively low molecular volume of cement, comp. 
Table I. However, in a finite heterogeneous system sol- 
vent flows tangential with respect to the interfaces are 
also expected and the argument given above is limited. 
In principle one has to consider the equations of hy- 
drodynamics for this purpose. However, this moving 
boundary problem is complex and we will not con- 
sider it in this form here. Instead we will determine 
the spatial and t emp oral solvent distribution v} aq ^ (t) 
according to Eq.(2.1). 



III. RESULTS 

From a general point of view cement hydration can 
be regarded as a heterogeneous (non equilibrium) solid 
phase transformation forming from an anhydrous solid 
phase two hydrated solid phases In contrast to 

known solid-solid phase transformations for example 
in metallic alloys, cement hydration bases strongly on 
the presence of a solvent phase (mostly water). The 
solvent controlls the transformation in a twofold sense: 
it is directly part of the chemical reactions and on the 
other hand it controlls the ion transport to a very 
large extent. 

The complexity of the problem and the features of 
the present approach are most convenient illustrated 
by Fig. [j]. Therein we show the solid volume frac- 



tions of anhydrous cement, e 



(C 3 S) 



hydrated cement, s 



(CSH) 



(left side), and of 



(right side), for three dif- 



ferent hydration times (a) t = 1500 s, (b) t — 2500 s 
and (c) t = 5000 s. The occupied volume fractions 
are color coded ranging from blue (0-20%), cyan (20- 
40%), green (40-60%) over yellow-orange (60-80%) to 
red (80-100%). The calculated micrographs Fig. [l]a- 
c show, as expected, spatially inhomogeneous nucle- 
ation and growth of CSH hydrate. The hydrate pre- 
cipitation and the associated growth of hydrate sur- 
face layers surrounding dissolving cement particles is 
relatively slow. However, in Fig. [j]b one especially ob- 
serves that the rate of CSH precipitation is spatially 
strongly varying (notice the red and cyan spots in the 
right image). This reflects the spatial fluctuations in 
chemical reactivities (local C3S surface/volume ratio), 
i.e., the reactivity of a 'fjord' is higher than that of a 
big 'lake'. Supported by experimental observations it 
has been argued that cement particle surfaces close 
to each other may indeed act as very strong inhomo- 
geneities leading to localized CSH nucleation forming 
'bridges' between adjacent particles Our micro- 
graphs confirm this picture. We note, that the pro- 
posed model does not include any (auto-)catalyticef- 
fect of CSH on its growth yet, compare Eq. (2.6). 
This point is intensively studied at the moment. 

As the hydrate precipitation/dissolution reaction 
rates depend on the local supersaturation values, the 
interesting question arises how global characteristic 
points of the cement hydration kinetics can be defined. 
In a recent experimental work it has been proposed to 
characterize different kinetic regimes occuring during 
the C3S hydration by means of parametric plots of 
'typical' calcium versus silicate concentrations as es- 
timated from calorimetric or conductimetric measure- 
ments ('kinetic path approach') ||. The experiments 
have been conducted employing stirred diluted C3S 
suspensions having a water/cement weight ratio be- 
tween 10 and 50. However, to obtain typical quanti- 
ties of such diluted suspensions employing numerical 
simulations one would have to consider very large sys- 
tems. Instead we will compare the experimental sus- 
pension data with numerical cement paste data having 
a water/cement ratio close to the practically impor- 
tant value 0.5. 

For convenience we reproduce in Fig. ^| some of the 
experimental data (o) of Ref. [)[ Initially one finds an 
'inductive' cement dissolution period (period AB in 



Fig. ||) . During this period the bulk solution is every- 
where undersaturated with respect to CSH and CH, 
the concentration gradients close to the cement parti- 
cles are very high, and ions rapidly distribute into the 
bulk solvent. We have observed in our simulations 
that CSH is already formed during the 'induction' 
period, however, to a very low extent. Because of 
the initially very high concentration gradients the ion 
products exceed the CSH solubility only in a very 
thin layer surrounding the cement particles. As the 
ions cummulate in solution the thickness of this pre- 
cipitating layer increases due to the broadening of the 
ion distributions. 

The precipitation counteracts the further increase 
in silica concentration both by chemical reaction and 
by decreasing the cement particles surface permeabil- 
ity. In result the silica concentration passes through 
a maximum (point B in Fig. ||), which might be in- 
terpreted as a global measure for the onset of CSH 
precipitation. During the period BC the silica con- 
centration decreases, while the calcium and hydroxyl 
concentrations continue to increase as the CSH pre- 
cipitation consumes only a fraction of these ions, see 
Eqs. (U) and (|1|. 

Typically after a couple of hours the solution be- 
comes oversaturated with respect to calcium hydrox- 
ide CH (point C in Fig. g). Both hydrates precipitate 
very slowly (experimentally several weeks between 
point C and D), lowering both calcium and hydroxyl 
concentrations. Point D corresponds to the state after 
infinite hydration time terminating the shown curve. 
The progress in time is indicated by arrows. Our nu- 
merical results (+) are in reasonable agreement with 
the experimental data (o). We have observed that 
calculated positions and values for point C depend 
a) on the employed initial water/cement weight ra- 
tio and b) on the three reaction rate constants k\ , &2 
and &3 in a rather complex manner. Contrary to this 
point B is found to depend only slightly on the rel- 
ative reaction rate constant k\jki over four orders of 
magnitude, see Fig. ga. A consequence of the much 
smaller solubility of CSH as compared to the second 
hydrate CH is that in the period AB only two chem- 
ical reaction are operative. Reaction Eq. (2.7) is in- 
operative. For k%/k2 ranging between 0.1 and 10 the 
maximum silica concentrations take constant values 
of about lmmolf liter while for k\/k2 — 10 3 we find 
« 3 mmol I liter . We have also performed simulations 
for similar kxjk-i, values employing different absolute 
reaction rate constants, see Fig. ga. We find a rela- 
tively good data collapse over four orders of magni- 
tude. This result is important insofar as the reaction 
rate constants have not been experimentally measured 
yet. In Fig. ||b we plot the momentaneous calcium 
concentration versus the silica concentration at point 
B for various reaction rate constants. The data are 
centered around a straight line of slope 1/3. Compar- 
ison of this result with the stoichiometric coefficients 



in Eq. ( |2.3| ) shows that the silica/calcium ratio of ions 
is determined by the dissolution reaction. 

The characterisation of the temporal advancement 
of the hydration process is an important problem. We 
show in Fig. |] some preliminary results for the time 
dependence of the mean volume fraction (e^ CSH ') of 



the cementious hydrate CSH. For 'moderate' reac- 
tion constants (o) the hydration is found to be rela- 
tively slow, i.e., after 70 hours of hydration we find 
only 10% CSH. This value is approximately two 
times lower than the experimental one observed by 
NMR 0. We note that both the hydration curve 
(o) in Fig. U and the parametric curve (+) in Fig. || 
belong to the same simulation. Apparently the hy- 
dration advancement for 'moderate' reaction rates is 
nearly constant in time, see curve (o) in Fig. ^. For 
comparison we also show in Fig. || a hydration curve 
for 'fast' hydration (+). After a few minutes of 'in- 
duction' period the hydration rapidly accelerates and 
goes already after 1 hour to completion (remaining 
C3S less than 1 %). This process is much to fast 
and it leads to typical ion concentrations of mol / liter 
which are unrealistic high. 

As long as the reaction rate constants have not been 
estimated from experiments yet one main difficulty in 
the modellization of cement hydration consists in find- 
ing appropriate values for the rate constants, which 
do not contradict experimental measurements of mean 
ion concentrations and hydration advancement. How- 
ever, the experimental results belong to hydration in 
space while our calculations correspond to two dimen- 
sional hydration. This point is being currently inves- 
tigated. 



IV. CONCLUSIONS 

We have presented a general, heterogeneous 
reaction-diffusion model for solid phase transforma- 
tion due to chemical dissolution/precipitation reac- 
tions. The model has focussed on the important in- 
dustrial problem of Portland cement hydration though 
it is more generally applicable to water/mineral sys- 
tems. 

We have tried to develop an 'open' approach based 
on physical and chemical considerations (mainly the 
laws of mass conservation and mass action). The 
model includes in its present form on a coarse grained 
length and time scale the full spatial distribution of 
solid and liquid phases, the three main chemical dis- 
solution/precipitation reactions, and the transport of 
ions due to diffusion. The proposed approach natu- 
rally allows to include for a varity of more or less im- 
portant phenomena depending on imposed conditions 
such as solvent flow effects, electrochemical effects, 
exothermic effects including heat conduction, pressure 
effects, inert filler effects etc.. 

We have tried to incorporate a reasonable amount 
of specific information about the cement hydration 
into the investigated model such as the stoichiome- 
try and kinetics of Portland cement dissolution, the 
precipitation/dissolution of the main cementious hy- 
drate (CSH) and of Portlandite (CH), their approxi- 
mate solubilities and molecular volumes, approximate 
values for the ion diffusivities in aqueous solution, and 
initial conditions for water immersed cement particles 
(sizes and spatial positions). 

The presented results demonstrate considerable 
richness and complexity of the cement hydration phe- 
nomenon close to controlled experimental situations. 



We have presented some calculated cement micro 
structures as they evolve in time. Nucleation and 
growth of hydrates is found to be strongly heteroge- 
neous in agreement with experimental observations. 
The problem of 'autocatalytic' effects on the precipi- 
tation processes needs further investigation. 

The presented parametric plots for the average con- 
centrations of ions in solution are in qualitative agree- 
ment with recent indirect experimental observations 
made for stirred cement suspensions §. The calcu- 
lated maximum silica concentrations and their cor- 
responding calcium concentrations are in reasonable 
quantitative agreement with experimental values. For 
high dissolution and low precipitation rate constants 
we find unacceptable high ion concentrations of or- 
der mol /liter. One way to circumvent this problem 
would b e th e introduction of a finite solubility for C3S 
in Eq. (^4|), however, such a solubility constant has 
not been experimentally determined yet. 

Furthermore we have investigated the variations of 
the maximum silica concentrations for various reac- 
tion rate constants. The concentrations are found to 
vary only slightly with the relative reaction rate con- 
stant fci/fe over four orders of magnitude. It would 
be very interesting to see how the second character- 
istic hydration point (point C in Fig. |^) depend on 
the initial water/cement ratio and on the employed 
reaction rate constants. 

Finally we presented first results on overall hydra- 
tion curves. This curves are not very realistic yet. It 
would be very helpfull to have experimental order of 
magnitude estimates for the three (unknown) reaction 
rate constants. 

For future work we are planning to conduct cal- 
culations in three dimensions with improved kinetic 
equations. 
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TABLE I. Stoichiometric numbers v\ 



(«) „(«) 



and 1/3 of the chemical reactions Eqs.(2.3), (2.5) and (2.7) respectively. 



The last three colums give the molecular volumes, the diffusion constants and the solubilities at room temperature and 
normal pressure. 



Species 
a 


la) 


Stoichiometric numbers 

(a) 


la) 

H 


Molecular volume 
[10~ 3 liter /mol] 


Diffusion constant 
[10 m r/s] 


Solubility 
[10~ 3 mol/liter] 


H 2 


-3 


-1 





18.01 







C 3 S 


-1 








73.1 (b) 





00 (fl 


C1.5SH2.5 





1 





80.40 (c) 





2.4 (d) 


CH 








1 


33.07 (a) 





20.2 (a) 


Ca 2 + 


3 


-3/2 


-1 




7.9 (e) 




OH- 


4 


-1 


-2 




53.0 (e) 




H 2 SiO\- 


1 


-1 







5.0 W 





(a) H.F.W. Taylor, Cement Chemistry, (Academic Press, London, 1990), pp. 125; {h) ibid, p. 15; (c) ibid, p. 152, 
we assume for the calculations a density of 2.35 g cm~ 3 ; ^ We have calculated this value from the supersaturation 
concentrations of calcium and silicate ions, as reported in Ref. |^; ' e ' R.A. Alberty and R.J. Silbey, Physical Chemistry, 
(Wiley&Sons, New York, 1992), pp. 839; ^' Measured values are not known to us. 




FIG. 1. The above figures show numerically calculated cement micro structures after (a) t = 1500 s, (b) t — 2500 s and 
(c) t = 5000 s. Images on the left side show the unhydrated cement phase (C3S) while those on the right side represent 
the hydrated, precipitated cement (CSH). The colors indicate the local solid volume fractions, i.e., blue (0-20%), cyan 
(20-40%), green (40-60%), yellow-orange (60-80%) and red (80-100%). The initial water/cement ratio is about 0.4, the 
reaction rate constants are ki = 10 -3 , ki = 10 -12 and k-j, = 10 -10 (in units of molmT 2 s _1 ). The employed cell size is 
Aa; = 10 -4 m, the linear system size L — lOOAa;, and the linear cement particle size £ = 5Ax. The time step is 0.01s. 
All other parameters as in Table I. 




(Ca 2+ ) [1(T 3 mol 

FIG. 2. Parametric plot of the evolution of silica versus calcium concentration ('kinetic path approach'); (o) exper- 
imental results after Ref. (+) numerical results. The employed reaction rate constants are ki = 5 ■ 10 -5 , = 10 -6 
and fe3 = 10 -8 (in units of molm" 2 s _1 ). The linear system size is L = 50Aa;, the cell size Aa; = 10 -4 m and the 
linear cement particle size I = 5Ax. The employed time step is 0.1s. Initial water/cement weight ratio 0.6. All other 
parameters as in Table I. 
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FIG. 3. Calculated silica ion concentration (H2Si0 2 ~) at the onset of CSH precipitation (corresponding to point 
B in Fig. ^|): (a) double logarithmic plot versus the relative reaction rate constant fei/fe, (b) plot against the mo- 
mentaneous (Ca 2+ ) ion concentration. Symbols: (o) ki = 10 -4 mol m~ 2 s" 1 , fei/fe ranging between 1 and 1000; (+) 
k\ = 10 -5 mol m~ 2 s -1 , fei/fe ranging between 0.125 and 500; (□) fci = 10 -6 mol m~ 2 s" 1 , fci/fe ranging between 0.167 
and 600; (x) fej = 10" 7 mo« m" 2 s" 1 , fci/fc 2 ranging between 0.167 and 0.5. Each point corresponds to an average over 
10 configurations having an initial water/cement weight ratio of about 0.5. All other parameters as in Fig. H. 
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FIG. 4. Double logarithmic plot of the average CSH volume fraction { e {CSH) ) versus hydration time t for two sets 
of reaction rate constants: (o) fci = 5- 10" 5 , k 2 = 10" 6 and fc 3 = 10" 8 ; (+) ki = 10" 2 , k 2 = 10" 16 and fc 3 = 10" 13 ; units 
in molm~ 2 s _1 . All other parameters as in Fig. @. 



